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Cosmic ray spectral hardening due to dispersion of source injection spectra 
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Abstract: The cosmic ray (CR) energy spectra measured with ATIC, CREAM and PAMELA showed that there is re- 
markable hardening for rigidity of several hundred GV. We propose that this hardening is due to the superposition of 
spectra from a population of sources, e.g., supernova remnants (SNRs), whose injection spectral indices have a disper- 
sion. Adopting proper model parameters the observational data can be well explained. It is interesting that the injection 
source parameters are similar with that derived from gamma-ray observations of SNRs, which may support the SNR- 
origin of CRs. Furthermore this mechanism provides an alternative explanation of the "ankle-cutoff" structure of the 
ultra high energy CR spectra. 
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1 Introduction 

The origins, acceleration processes, propagation and inter- 
actions of cosmic rays (CRs) are still open questions, even 
after about one century since the discovery of CRs. It is 
generally believed that Galactic CRs (GCRs) are acceler- 
ated by the astrophysical shocks like supernova remnants. 
After the production, GCRs are then injected into the in- 
terstellar space and propagate diffusively in the Galactic 
magnetic field. During the propagation, CRs will interact 
with the interstellar medium, radiation field and magnetic 
field, which result in spallation and energy loss of the parti- 
cles, and the production of secondary particles (see [ 1 1 for 
a review). 

In recent years more and more accurate data of the CR 
spectra and composition are available. The balloon-borne 
experiment Advanced Thin Ionization Calorimeter (ATIC) 
measured the CR spectra of various species of nuclei and 
showed the deviation from single power-laws of the spec- 
tra fl2] [3) . ATIC data also revealed the difference between 
proton and Helium spectra. Another balloon experiment 
Cosmic Ray Energetics And Mass (CREAM) measured the 
energy spectra of the major species from proton to iron 
with relatively high precision, and reported a remarkable 
hardening of the spectra of most heavy species at ~ 200 
GeV/nucleon HE]. Such spectral features were confirmed 
most recently by the satellite experiment Payload for Anti- 
matter Matter Exploration and Light-nuclei Astrophysics 
(PAMELA). PAMELA measured the proton and Helium 
spectra up to rigidity 1.2 TV with high precision, and the 



data show clearly a spectral hardening at rigidity ~ 200 
GV 0, which is basically consistent with the results of 
CREAM and ATIC. The spectral indices of proton and He- 
lium were also found significantly different. 

The hardening of the CR spectra challenges the conven- 
tional idea that the spectra below the so-called "knee" are 
simple power-laws. More complicated scenarios of CR ori- 
gin, acceleration or propagation are needed to explain the 
data. Models possibly to explain such a spectral hardening 
include the multi-component sources [7|, or the nonlinear 
particle acceleration scenarios where the feedback of CRs 
on the shock is essential (e.g., If8l l9l ITOl ITT1 ^ . 

In Ifl2l we propose that the hardening of the observed CR 
spectra is due to the dispersion of the injection spectra of 
a source population such as supernova remnants (SNRs). 
The basic fact that the superposition of a series of spectra 
with different power-law indices would result in an asymp- 
totic hardening of the final spectra of CRs was recognized 
long ago [13). Such an effect was also employed to explain 
the "GeV excess" of Galactic diffuse 7-rays observed by 
EGRET lfl4lfT5l. 

The observations of various kinds of candidate high en- 
ergy CR sources such as SNRs, active galactic nuclei and 
gamma-ray bursts, show indeed there is significant dis- 
persion of the source spectra. For example, the spectral 
modeling of the 7-ray emission from several SNRs ob- 
served by Fermi and ground-based Cherenkov telescopes 
show that the low and high energy spectral indices are 
71 « 2.15 ± 0.33 and 72 2.54 ± 0.44, assuming a 
hadronic scenario of the 7-ray emission 1 16 1. In the follow- 
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ing we will show how such a scenario can natually explain 
the observational data of the CR spectra. 

2 Model and result 

We assume that GCRs are originated from SNR-like 
sources. The injection spectrum of each source is assumed 
to be a broken power-law function of rigidity. The spec- 
tral indices are assumed to be Gaussian distributed around 
the average values. The break rigidity Rf, is shown to be 
several to tens of GV fl6l . Here we assume the logarithm 
of Bb is randomly distributed in some range. The normal- 
ization of each source is derived assuming the same total 
energy of CRs above 1 GeV for all sources. The detailed 
parameters are given below, according to the fit to the ob- 
servational CR data. 

The CR propagation is calculated with the public GAL- 
code 111711 . In this work we adopt the diffusive reac- 
celeration frame of the propagation model, with the main 
parameters D = 5.8 x 10 28 cm 2 s -1 , 6 = 0.33, v A = 32 
km s" 1 and Zh — 4 kpc. The input source spectra are 
adopted as the superposed ones of many sources. The B/C 
ratio is found to be well consistent with data and insensitive 
to the source spectrum. 

For the low energy spectra (below ~ 30 GV), the force- 
field approximation is adopted to model the solar modula- 
tion effect [18|. The modulation potential $ varies from 
— 200 MV to - 1400 MV, depending on the solar activity. 
The modulation potential was estimated to be 450 — 550 
MV for the time when the proton and Helium data used 
in this work were recorded by PAMELA (6) . In this work 
we adopt $ = 550 MV for proton and Helium. For Car- 
bon, Oxygen and Iron nuclei the low energy data used are 
from HEA03, for which we adopt $ = 750 MV. A higher 
modulation potential for HEA03 was also found in |[T9l . 

We further adopt a break/cutoff at high rigidities (^PV) to 
approach the "knee" of the CR spectra. The physical mod- 
els of the knee include the propagation/leakage effect from 
the Galaxy, or interactions with background particles [20|. 
Phenomenologically we adopt two kinds of cutoff/break 
to model the knee structure of the total spectra: a sub- 
exponential cutoff case with the energy spectrum above 
the injection break i?~ 72 exp [— (i?/i? c ) 1//2 ], and a broken 
power-law case with energy spectrum above the injection 
break Rr l2 {l + R/R^)~ l , In both cases we assume the 
break rigidities R c of different nuclei are the same, i.e., 
the break energies are Z-dependent. The parameter R c is 
adopted to be a constant instead of that varying for differ- 
ent sources. This assumption is reasonable for the propa- 
gation/leakage and interaction models. However, we may 
note that for the acceleration limit models, the break en- 
ergies of different sources should have a dispersion. It is 
tested that the results with dispersion of the break rigidities 
are similar to that with a proper constant R c . 

The calculated energy spectra of proton, Helium, Carbon, 
Oxygen, Iron and the total spectrum, together with the ob- 



servational data are shown in Figure Q] The parameters 
used in the calculation are compiled in Table Q] We find 
relatively good agreement between the model calculation 
and the observational data. It is interesting that the injec- 
tion source parameters are similar to those inferred from 
the 7-ray observations of SNRs. This might be evidence 
that SNRs are the sources of GCRs below ^PV. 

Note that our model prediction is a gradual hardening of 
the CR energy spectra, which seems to be consistent with 
the overall structures of data in a wide energy range. How- 
ever, if we focus on the detailed structures revealed by in- 
dividual experiment, we may still find some inconsistency. 
For example the PAMELA data show a very sharp dip at 
~ 200GV and a gradual softening below the break rigid- 
ity J6). The CREAM data indicate that all species of nu- 
clei experience a hardening at ~ 200 GeV/n [5 |. The 
very fast break of the spectra shown in both PAMELA 
and CREAM data cannot be well reproduced in the present 
model. Before giving a conclusive judgement about this 
issue, we may need future better measurements of wide en- 
ergy band spectra by e.g., the Alpha Magnetic Spectrom- 
eter (AMSO20) and the Large High Altitude Air Shower 
Observatory (LHAASO, ED). 

3 Impact on secondary particles 

The hardening of the primary CR particles should have im- 
prints on the secondary particles such as positrons [22], 
diffuse 7-ray s, and antiprotons 11231 . In ll24l a systematic 
study of the secondary particles including B/C, antipro- 
tons and diffuse 7-rays in different scenarios of the spectral 
hardening was performed. As an example, here we calcu- 
late the predicted hadronic-origin diffuse 7-ray fluxes in 
this scenario. The 7-ray production spectrum from pp in- 
elastic interactions is calculated using the parameterization 
given in 11251 . The contribution to 7-rays from heavy nuclei 
in both the projectile and target particles is approximated 
with a nuclear enhancement factor cm = 1.84 [26 1 . We 
show in Figure [2] the ratios of the hadronic component of 
the diffuse 7-ray fluxes between our model expectation and 
that of the traditional single power-law CR spectrum. It is 
shown that the 7-ray flux will also experience a hardening 
above ~ 50 GeV. The most remarkable hardening effect 
will be in the very high energy range (>TeV), which is out 
of the capability of Fermi telescope. 

Since the total diffuse 7-rays consist with hadronic, lep- 
tonic and the extra-galactic components, the hadronic one 
is not a direct observable. As we know that in the Galac- 
tic plane the diffuse 7-rays should be dominated by the 
hadronic component [27], we expect such a predication can 
be tested with high precision diffuse Te V 7-ray observation 
of the Galactic plane. 

For positrons and antiprotons we will expect similar re- 
sults, however, the quantitative effects should depend on 

1. http://galprop.stanford.edu/ 

2. http://ams.cern.ch/ 
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Figure 1: Energy spectra of proton (top-left), Helium (top-middle), Carbon (top-right), Oxygen (bottom-left), Iron 
(bottom-middle) and the all-particle one (bottom-right). The solid line in each panel represents a sub-exponential cut- 
off behavior of the high energy spectra around the knee region, while the dashed line is for broken power-law type. For 
references of the data please refer to fl2l . 



Table 1: Source parameters: injection spectra 71, 72 and break rigidity high energy cutoff rigidity R c and solar 
modulation potential $. 
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the source distribution and propagation model. The current 
antiproton data from PAMELA can not probe such a hard- 
ening effect yet. Also the PAMELA positron excess above 
several GeV should not be originated from the hardening 
of the CR proton and Helium spectra. 

4 Alternative explanation of the ultra high 
energy CR spectra? 

We also note that the dip-cutoff structure shown in Figure 
[T]is very similar to the ankle and Greisen-Zatsepin-Kuzmin 
(GZK) cutoff structure of the ultra high energy CRs (UHE- 
CRs). If the UHECRs are originated from a population 
of sources instead of a major one, it is expected that the 
UHECR spectra should also suffer from such a hardening 
due to the superposition effect. It might provide us an alter- 
native point of view to understand the observational energy 
spectra of UHECRs. 

We show an illustration of the predicted UHECR spectra, 
together with the HiRes data [28], in Figure [3] Here the 
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Figure 2: Ratio of the hadronic component of the diffuse 
7-rays between the dispersion scenario expectation and the 
single power-law model. 

chemical composition of UHECRs is assumed to be pure 
protons. The injection spectrum of UHECRs is adopted to 
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be a broken power-law function with an exponential cut- 
off 1291 . The logarithm of break energy log(.Ef,/eV) is as- 
sumed to be uniformly distributed in [17, 18], and the spec- 
tral index is 2.0 ± 0.2 below E b and 3.6 ± 0.6 above E b . 
The cutoff energy E c is adopted to be 5 x 10 19 eV. The 
result shows a good description of the observational data. 

If the UHECRs are produced at cosmological distances, 
however, the interactions between UHECRs and the cos- 
mic background photons are unavoidable l29l . Therefore 
if the superposition effect is responsible for the shape of 
the UHECR spectrum, we may in turn expect that UHE- 
CRs are produced locally or even in the Galaxy ll30l . More 
generally, it is possible that both the superposition and the 
interaction effects are in operation to give the ankle-GZK 
structure of UHECRs. 




E (GeV) 



Figure 3: Calculated energy spectra of UHECRs for pure 
protons, compared with the HiRes data |28|. 



5 Summary 

In summary we propose that the superposition of the en- 
ergy spectra from CR sources with a dispersion of injection 
spectra is responsible for the recently reported hardening 
of the CR spectra. If the CRs are indeed originated from 
a population of sources instead of a single major source 
(e.g., iFJTI ). such an asymptotic hardening effect due to dis- 
persion of source properties is inevitable. It is interesting to 
note that the source parameters derived are similar to that 
inferred from 7-ray observations of SNRs, which might 
support the SNR-origin of GCRs below the knee. Finally 
we discuss the same mechanism as an alternative explana- 
tion of the ankle-GZK structure of UHECR spectra. 
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